Stress and large-scale spatial structures in dense, driven granular flows 
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We study the appearance of large-scale dynamical heterogeneities in a simplified model of a driven, 
dissipative granular system. Simulations of steady-state gravity-driven flows of inelastically colliding 
hard disks show the formation of large-scale linear structures of particles with a high collision 
frequency. These chains can be shown to carry much of the collisional stress in the system due to a 
dynamical correlation that develops between the momentum transfer and time between collisions in 
these "frequently-colliding" particles. The lifetime of these dynamical stress heterogeneities is seen 
to grow as the flow velocity decreases towards jamming, leading to slowly decaying stress correlations 
reminiscent of the slow dynamics observed in supercooled liquids. 



PACS numbers: 81.05.Rm, 45.70.-n, 83.10.Pp 



I. INTRODUCTION 

Granular materials are ubiquitous in nature and play 
an important role in many industrial applications, yet 
the appropriate coarse-grained theories for the descrip- 
tion of these materials remain elusive due to the rich phe- 
nomenology they display in response to external pertur- 
bation P, 12 • A common example of the unusual nature 
of granular systems is the phenomenon of force chains, 
highly-stressed extended linear clusters first seen in ex- 
periments on static granular piles. These striking 
inhomogeneities are thought to be responsible for a host 
of observed nonlinear effects in static systems including 
interesting features in distributions of contact forces be- 
tween particles Q. However, the reason for the existence 
of the force chains has not been determined; one possible 
explanation is that correlations develop while the mate- 
rial is flowing which then get " frozen in" when the system 
comes to rest. 

Spatial structures have been directly observed in ex- 
periments on flowing granular materials More re- 
cently, measurements of velocity correlations in the sur- 
face layer of particles flowing down an inclined plane have 
yielded a length scale which appears to grow as the flow 
is arrested [6j. While detection of such dynamical het- 
erogeneities in numerical studies has been more difficult, 
there is some evidence for their existence in simulations 
of gravity-driven flow through a hopper and in im- 
ages of the contact force network in simulations of chute 
flow |8( ■ The presence of such structures will play a cru- 
cial role in determining the dynamics of the flowing sys- 
tem, especially as the transition to a static system is 
approached. The development of spatial and temporal 
correlations near the transition to the static state must 
be investigated carefully in order to understand the na- 
ture of the transition from the flowing liquid state. 

Dynamical heterogeneities have been observed in other 
systems which exhibit a similar slow dynamics as the 
transition from fluid to solid is approached, such as col- 
loids near the glass transition. Fast-moving particles 
were observed to be spatially correlated with a charac- 
teristic cluster size that increased as the glass transition 



was approached :9j. In supercooled liquids, spatial inho- 
mogeneities can be identified via a time-dependent four- 
point (two-time, two-space) correlation function |10| . 
The interesting commonalities between these disparate 
systems have been explored in the context of the jam- 
ming phase diagram proposed in Refs. Il2j . ft should 
be noted that in these studies the focus has been primar- 
ily on the search for a common static signature (i.e. the 
distribution of contact forces P(f)). Thus the proper 
framework to describe the dynamics of systems which 
have a jamming transition remains an open question. 
If the presence of dynamical heterogeneities governs the 
dynamics of these far-from-equilibrium systems near the 
jamming transition, then while the microscopic process 
by which these structures form may vary between differ- 
ent types of systems, some unified dynamical description 
may be applicable once the structures have formed. The 
"soft glassy rheology" (SGR) model [H 03 describes 
soft, nonergodic systems such as foams and dense emul- 
sions and has successfully predicted many of the exper- 
imentally observed results. Application of this type of 
"trap- like dynamics" to granular systems has also had 
some promising initial results |l5j . 

In this paper we explore the connections between dy- 
namical heterogeneities and slow dynamics in a simpli- 
fied model which allows us to focus on the essential ef- 
fects of driving and dissipation in a flowing granular sys- 
tem, while still reproducing observable results from re- 
lated experimental systems [l6| . Simulations have been 
performed of a two-dimensional gravity driven system of 
frictionless, bidisperse hard disks in a hopper geometry. 
The disks undergo instantaneous, inelastic binary colli- 
sions and propagate under gravity in between collision 
events. Despite the simplified dynamics of the simu- 
lation, several interesting collective effects arc seen to 
emerge. As in previous studies 0, extended linear struc- 
tures defined by inhomogeneities in the spatial distribu- 
tion of the collision frequency are visible, and their pres- 
ence has a measurable effect on distributions of impulse 
and time intervals between collisions for a given parti- 
cle. In this study we observe that the particles which 
make up the linear structures display a dynamic corre- 
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lation between their impulse and collision time which is 
reflected in measurements of the collisional stress tensor; 
particles identified as part of the structures appear to ex- 
perience higher than average collisional stress compared 
with the remainder of the system. Using this measure of 
collisional stress we can probe further into the nature of 
these structures and their relevance to the dynamics of 
the system. 



II. SIMULATIONS 

The grain dynamics used in the simulations are sim- 
ilar to Ref. • Specifically, (i) at each interparticle 
collision, momentum is conserved but the collisions are 
inelastic so the relative normal velocity is reduced by the 
coefficient of restitution p,; (ii) to allow the side walls to 
absorb some vertical momentum we impose the condition 
that collisions with the walls are inelastic with a coeffi- 
cient of restitution fx wa ii in the direction tangential to the 
wall and (iii) since we wish to observe the system over 
many events, particles exiting the system at the bottom 
must be replaced at the top to create uniform, sustained 
flow. Note that collisions are instantaneous and in be- 
tween collisions the particles are driven by gravity. To 
avoid the phenomena of inelastic collapse the coefficients 
of restitution /i and fJ, wa ii are velocity dependent; if the 
relative normal velocity between particles or between a 
given particle and the wall is less than some cutoff v cu t 
then the collision is presumed to be elastic [l^- The 
flow velocity is controlled by adjusting the width of the 
hopper opening. We also introduce a probability of re- 
flection p at the bottom which reduces the time needed to 
reach steady-state flow. Typically, our simulations were 
done on bidisperse systems (diameter ratio 1:1.2) of 1000 
disks, with /i = 0.8, \i wa ii = 0.5, v cu t = 1 x 10~ 3 and p = 
0.5. The simulation was run for a total time of 1 x 10 3 in 
simulation time units (smaller particle diameter d s and 
gravitational constant g are both set to 1) with the initial 
time interval of 5 x 10 2 discarded before recording data to 
ensure the system has reached steady state. During the 
total time interval of 500 over which we are evaluating 
the data, a given particle will pass through the hopper 
5-10 times depending on the flow velocity. 

As the width of the hopper opening is decreased the 
average flow velocity decreases, and at some minimum 
opening sustained flow is no longer observed. Close to 
this width the mass flow exiting the hopper appears to 
become very intermittent, with large outflow of mass oc- 
curing on short timescales followed by long time intervals 
where few particles exit the system. Intermittent dynam- 
ics has also been observed in simulations of flow down an 
inclined plane where the angle of incline is approximately 
equal to the angle of repose of the grains [8( . Addition- 
ally, some preliminary measurements of spatial velocity 
correlations in our hopper simulations parallel to the flow 
direction indicate long range correlations on length scales 
comparable to the system size for slow flow velocities . 



Spatial velocity correlations with long length scales are 
also seen on the surface of incline plane flow [(| . 



III. COLLISION FREQUENCY AND 
LARGE-SCALE STRUCTURES 





FIG. 1: Images of the system with particles coloured ac- 
cording to the collision frequency in the last time interval 
At = 0.02 for v f = 0.86 (left) and v f = 2.2 (right), (v) is 
defined as the mean collision frequency (average taken over 
all particles) for that time interval. 

In studies of two-dimensional freely-cooling gases [2(j 
some evidence of dynamical large-scale structures form- 
ing in the system was first observed by identifying all 
particles experiencing inelastic collapse as the assembly 
of grains cooled. These particles formed linear structures 
within dense clusters and experienced collision frequen- 
cies significantly higher than that of the other particles 
in the system. Due to the effect of driving in our flowing 
system, we do not observe inelastic collapse. However, 
we can ask an analogous question: "How many collisions 
does a given grain undergo in a given time interval At?" 
The collision frequency for a particular particle i in At is 
given by Vi — Nij At where N\ is the number of collisions 
experienced by particle i in At. We can construct images 
of our simulated system at regular intervals and colour 
individual disks according to the value of the collision fre- 
quency for that disk. As can be seen in Fig. 2] particles 
with high collision frequencies form linear chains reminis- 
cent of the structures observed in 2D freely-cooling gases 
and the transient solid chains postulated by the hydro- 
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dynamic model of Ref. 21] . These chains often appear 
to terminate at the boundary of the system, forming an 
arch-like structure in the centre of the hopper. From the 
pictures it is also evident that collision direction is typ- 
ically oriented along the chain direction (i.e. particles 
within a chain collide with other particles in the chain 
and not with neighbouring particles which arc not part of 
the structure) . Thus the velocities of frequently-colliding 
particles are highly correlated in direction, and we can 
evaluate the stability of this structure with respect to 
perturbations arising from collisions with other (uncor- 
related) particles. Consider an initially isotropically dis- 
tributed dense assembly of completely correlated disks 
and perturb one disk by giving it a velocity v. The col- 
lision rule for two inelastic hard disks i,j is [Fjj: 



(1) 







Vi 


1+fl 


-Vi • q + Vj • q 


q 
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v ; • q - Vj • q 





where q is a unit vector along the centre of mass line of 
the colliding disks and v' denotes the velocity of a given 
particle after the collision. Rewriting Eq. ^ in terms of 
the relative velocity of particle j to particle i before the 
collision Sv(t) = vj — Vi and after the collision Sv(t + 
At) = vj — vj, we can obtain an approximate expression 
for the time derivative of the relative velocity in the limit 
of At 0: 



dSv(t) 
dt 



At 



[-Jv(t) • q] q 



(2) 



Define qu and q± as the components of q parallel and 
perpendicular, respectively, to the initial perturbing ve- 
locity. From Eq.|21it is easy to see that (-3P-) — —C q (qn) 
and (^} ~ -C q {q ± ) K where C q = ^ [-*v(t) • q]. 
However, q± is equally likely to be positive or negative 
thus (q±) ~ while qu is always positive and so (q\\) > 0. 
Substituting this result into the above expressions for the 



/ dSv 



time derivatives yields (—^f-) < and ( ^ ) ~ indicat- 
ing that velocity differences parallel to the initial pertur- 
bation will decay while velocity differences perpendicular 
v will remain. Given this, and the organizing influence of 
the inelastic wall collisions leading to velocities becoming 
normal to the wall, the appearance of collisional chains 
anchored to the wall becomes plausible. 



IV. STRESS 

The unusual rhcology of granular matter in response to 
external stress, as well as the phenomena of force chains 
prompts us to investigate stress in this simple model of 
flowing granular matter. Specifically, are the frequently- 
colliding chains described in the previous section the dy- 
namical analog of force chains? We can explore this pos- 
sible connection by measuring the stress in our system 
and determining if a correspondence exists between the 





FIG. 2: Images of the system with principal axis correspond- 
ing to the maximum eigenvalue A m of the collisional stress 
tensor defined during the last time interval At = 0.02 for 
t)f = 0.86 (left) and Vf = 2.2 (right) plotted at the particle 
positions. The principal axis line is coloured and its width 
scaled according to the magnitude of the maximum eigen- 
value A m (the thicker the line the higher the value of A m ). 
Thin black lines indicate the principal axis of particles ex- 
periencing less than the minimum stress plotted (0.3(A) for 
these pictures). Note that these pictures were constructed for 
the same time interval as the images in Fig. Q (A) is defined 
as the mean A m . (average taken over all particles) for that 
time interval. 



collisional chains and the spatial distribution of stress. 
Additionally, we can potentially use the measured stress 
to calculate relevant time and length scales associated 
with the structures, as well as making a clearer con- 
nection between our model and other systems with soft, 
glassy rheology. 

As described in Ref. (2^ > macroscopic stress in discrete 
particle systems develops as a result of two microscopic 
mechanisms of momentum transfer: (a) transport of mo- 
mentum due to the fluctuations of the individual particle 
velocities (the kinetic stress au) and (b) transport of mo- 
mentum at interparticle collisions (the collisional stress 
cr c ). The total stress tensor is then the sum of the kinetic 
and collisional stresses. However, for the densities we ob- 
serve, we find as in Ref. [22| that the magnitude of the 
kinetic stress contribution is several orders of magnitude 
lower than the collisional stress. Thus we will consider 
only the collisional stress contribution in the analysis to 
follow. For a given particle i, the \xv component of the 
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collisional part of the stress tensor at time t is given by: 



where /z, are Cartesian coordinates. The sum is taken 
over all collisions a experienced by particle i during time 
t — At — > t, where the interval At is chosen such that 
(t) <C At <C 1/vf where r is the time between successive 
collisions for a particular particle. This constraint en- 
sures that the particle being evaluated experiences many 
collisions during At but is not significantly rearranging 
relative to its neighbours. There is a separation of time 
scales (which will be demonstrated more concretely in the 
analysis to follow) that ensures an interval At satisfying 
this constraint is easy to locate. Note that for hard-disk 
collisions I a will always be along r^, therefore the I Q ri Q 
term in the sum is simply (d,; + dj)I a /2 where di, dj are 
the diameters of particle i and the other particle involved 
in collision a respectively. 

Similarly to the images of the system at a given time t 
constructed for the collision frequency, we can make com- 
plementary images of a c (t) (Fig. |5J). For every particle 
we calculate a c {t) as described in Eq. [31 and extract the 
maximum eigenvalue X m along with the corresponding 
principal axis. A line along the direction of this principal 
axis is plotted at the particle position, and then coloured 
according to the value of A m relative to the average value 
of the maximum eigenvalue measured in that time inter- 
val (A). Fig. [2 shows some sample images constructed 
in this way for the same systems and times as in Fig. ^ 
Comparison of the two figures reveals a strong correla- 
tion between the frequently colliding particles and the 
highly-stressed particles. Additionally, the principal axis 
of the stress for the highly-stressed/frequently-colliding 
particles is typically aligned along the chain direction. 
There is also a lack of a obviously growing length scale 
as Vf — > 0, however, timescales associated with these 
highly-stressed structures show a significant change with 
flow velocity as discussed below. 

Note that the value of X m will depend on the ratio 
^2 I a /At, and thus, since the frequently colliding par- 
ticles experience many collisions in At (i.e. i/, « 1/t > 
(v), where r is the time between successive collisions 
for a particular particle), then the sum over I a can be 
large even if the average impulse per collision is small. 
Thus A TO (and similarly a c (t)) will depend crucially on 
the ratio of impulse / and collision time r, and any cor- 
relation between these quantities. Therefore, to gain a 
better understanding of the correspondence between the 
frequently-colliding particles and the highly-stressed par- 
ticles we have studied the joint distribution of impulse 
and collision time P(I,t). 



V. DISTRIBUTIONS OF IMPULSE AND 
COLLISION TIME 

In previous work |7j we modelled recent experiments 
on dense, gravity-driven monodisperse granular flows [l6|. 
Both the experiments of Ref. ^(| and the simulations we 
performed measured the distribution of impulses trans- 
ferred at each collision, P(I) as well as the distribution of 
time intervals between the instantaneous collisions P(t). 
In the experiment the measurement was made at the wall 
by using a transducer positioned a short distance up- 
stream from the hopper opening to detect the impulse 
transfer as a function of time to that point on the wall. 
In the simulation we can separately analyze data from 
events located within the bulk, at the walls or at a spe- 
cific point on the wall to mimic the experimental mea- 
surement. For most of our discussion we will focus on 
observations of distributions taken from events occurring 
in the bulk of the material. 

Both the experiment and the simulation revealed the 
same effects in the measured distributions: (i) an increase 
in the number of small impulse events as the flow veloc- 
ity is decreased as well as (ii) an increase in the number 
of small collision time events as the flow velocity is de- 
creased. In the impulse distribution in particular, this 
change could be directly linked to the appearance of the 
frequently-colliding clusters pictured above. Addition- 
ally, a simple one-dimensional model was devised that 
demonstrated how the existence of linear clusters of par- 
ticles which collide primarily along the chain direction 
would lead to the observed behaviour in P(I) 7J. The 
flow-rate invariant exponential form of P(I) at large im- 
pulse as measured in both experiment and simulation was 
shown to be a reflection of the local velocity distribution, 
which had an exponential tail (note that the impulse dis- 
tribution of uncorrelated particles is essentially the con- 
volution of the individual velocity distributions). Using 
the same ID toy model we were able to demonstrate that 
the large impulse region of P(I) is determined primarily 
by the form of the velocity distribution as observed. 

In the bidisperse system, the picture retains many of 
its original features with some added complexity. Defin- 
ing the scaled impulse / = //(/) where (I) = J IP(I)dI 
we can consider the behaviour of the joint distribution 
of impulse and collision time P(I,t). Shown in Fig. [2k- 
as a function of / for two values of r, several interesting 
features are seen to emerge. For both values of r, P(I, r) 
still shows an increase in height at impulses smaller than 
the average impulse as the flow velocity is decreased. 
However, for r = 10~ 3 the large impulse region of the 
distribution is no longer flow rate independent as in the 
monodisperse system, and an increase in height at im- 
pulses much larger than the average is also visible [23| . 
As before, this large impulse behaviour can be linked to 
an accompanying shape change in the velocity distribu- 
tion (Fig. I3t>), and is similar to the behaviour of P(J) in 
three-dimensional MD simulations of soft spheres [24| . 

If we look at P(I, r) as a function of / for a given flow 
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FIG. 3: (a) Joint impulse-collision time distribution P(I, r) 
at varying flow velocities. The inset shows the small im- 
pulse region of the distribution in more detail. Note that 
the two cuts of P(I, t) have been separately normalized for 
easier comparison, (b) Corresponding velocity distributions 

P(Vy/{Vy)). 



velocity, it is evident that the impulse and collision time 
are correlated (i.e. small impulses tend to be associated 
with small collision times). To clarify the roles played 
by the frequently-colliding and rarely-colliding particles 
in the observed shape changes in P(I, r) we can divide 
the joint impulse-collision time distribution into contri- 
butions from these two populations (shown in Fig. for 
Vf = 1.13). For this purpose, based on review of the 
images we chose a threshold frequency of 5(f) where 
(v) = (1/iV) J2i v i an d defined all particles with collision 
frequency Vi above this threshold as frequently-colliding. 
Rarely-colliding particles are similarly defined as parti- 
cles with < v i < 0.5(f). It appears that the corre- 
lation between impulse and collision time is stronger for 



FIG. 4: Impulse distribution P(I, r) showing the separate 
contribution from the frequently-colliding (FC) particles (f j > 
5(f)) and the rarely-colliding (RC) particles (yi < 0.5(f)) at 
v f = 1.13. 



the frequently colliding particles, and that these particles 
are the primary contributors to P(I, r) at small I and 
t [2j| . Additionally, as expected from the observations 
made on the large impulse region of the total P(I,t), 
at longer collision times and large impulses the shape of 
the distribution is governed by the contribution from the 
rarely colliding particles (which in turn reflects the shape 
changes in the velocity distribution shown in Fig. EJj) . 

The impulse/collision time correlation observed in 
P(I, t) is further evidence that velocities of frequently- 
colliding particles are correlated with small fluctuations 
leading to collisions with I < I. Note that this cor- 
relation in the velocities is apparent even in the ab- 
sence of any density inhomogeneity. The connection be- 
tween the frequently colliding particles and the highly- 
stressed particles is not surprising in the context of the 
impulse/collision time correlation discussed above, which 
would give rise to values of A m ^> (A). Therefore the 
stress heterogeneities seen in Fig. [21 are directly related 
to the dynamic correlation reflected in P{I,t). Given 
that this correlation is seen to become stronger as the 
flow velocity decreases, it is reasonable to expect some 
kind of flow rate dependent behaviour in time and length 
scales associated with the collisional stress. As previously 
described, there does not appear to be any significant 
change in length scale of the structures with flow veloc- 
ity. As a first attempt to quantify a trend in the time 
scale, we will calculate the autocorrelation C\(t) of A m 
and determine the dependence of the time scale for the 
decay of C\{t) onvf. 
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VI. RELAXATION OF STRESS 
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FIG. 5: Autocorrelation of maximum eigenvalue A m of a c (t). 
Inset shows the time scale ro.i at which (A m (f)A m (0)) = 0.1 
plotted against l/v/. 

The autocorrelation of the maximum eigenvalue 
C\(t) = (A m (i)A m (0)) of A m is calculated as: 



(Xm(t)Xm(0)} 



NN t0 



where (A m ) is the time- averaged value of the maxi- 
mum eigenvalue and 7V to is the number of time origins 
to. C\(t) (shown in Fig. [S] for varying flow velocities) ap- 
pears to have three decay regimes: (i) a short power-law 
regime at small time t ~ (r), (ii) a longer slow (possibly 
logarithmic) decay at intermediate time (r) <C t <C 1/vf 
and (iii) a more complicated decay at long time t > l/vf 
which involves a shallow negative dip. This three-stage 
form of C\(t) is reminiscent of density-density autocor- 
relations measured in supercooled liquids at low temper- 
atures |26| , which show a power law regime at short time 
scales, a plateau at intermediate time scales and a decay 
to zero at much longer time scales. This implies that the 
intermediate slow decay observed in the granular system 
may be resulting from a type of "caging effect" intro- 
duced by the correlated motion of particles participating 
in one of the highly-stressed/ frequently-colliding chains. 
This type of C\(t) indicates a growing separation of time 
scales as the flow velocity is decreased, which is an inter- 
esting similarity between granular flows and supercooled 
liquids. 

We can get a measure of the time scale associated with 
this intermediate slow decay of the stress by extracting 
the time to.i at which C\(t) has decayed to 10% of its 
original value (i.e Ca(tq.i) = 0.1). Plotting tq.i vs 1/vf 
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FIG. 6: (a) Autocorrelation of maximum eigenvalue A m of 
<j c (t) for both the frequently-colliding (FC) particles and the 
whole system at Vf = 2.05 and Vf = 0.916. Inset shows the 
time scale ro.i at which (A m (t)A m (0)) = 0.1 plotted against 
l/i)f for both the FC particles and the whole system, (b) 
Stress autocorrelation for ID weakly driven inelastic granular 
gas. 



(see inset to Fig. |SJ) clearly indicates that ro.i diverges 
as l/v f. If we associate this time scale with the life- 
time of the highly-stressed spatial structures visible in 
Fig. |21 then one can see that as the flow velocity de- 
creases the structures last for increasingly longer times. 
C\(t) for the frequently-colliding particles only is shown 
in Fig. [HJi. These curves were calculated by choosing only 
those time origins to where a given particle had been iden- 
tified as frequently-colliding according to the criteria de- 
fined above (z^ > 5(f)). Note the particle is not required 
to be frequently-colliding during the entire time interval 
over which C\(t) is calculated, thus we expect the long 
time behaviour of these curves to be similar to that of 
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the entire system. The same three-regime behaviour is 
evident, but the time scale of the decay r^f is longer 
(tq i = 2t .i). However, the scaling of the decay time 
with the flow velocity is the same (see inset to Fig. EJi). 

The role of the frequently-colliding/highly-stressed 
particles in determining the shape of the autocorrela- 
tion can be clarified by considering a related system; a 
one-dimensional assembly of particles which experience 
instantaneous, binary, inelastic collisions. To establish a 
steady-state, a fraction of the particles are weakly driven 
in between collision events. This ID system represents an 
idealized, isolated analog to the frequent-collision chains, 
and thus the stress autocorrelation measured in this sys- 
tem should be representative of the stress relaxation be- 
haviour of these structures (note that a c (t) in the ID 
system is a scalar quantity) . C\(t) for the ID system is 
shown in Fig.m> 27| . From this plot it can be seen that 
the initial stress relaxation is slow and possibly logarith- 
mic, similar to the slow decay regime of C\(t) in the 2D 
flowing system. From this result it appears that the ob- 
served dynamical heterogeneities do dominate the stress 
relaxation in the flowing granular system at intermediate 
time scales (t) C i < 

VII. CONCLUSIONS 

The picture which is emerging from these simulations 
is that large-scale highly-stressed structures analogous to 
force chains in static systems can form even in a sim- 
plified model of inelastic hard disks. These structures 
experience a slow relaxation from collisional stress at in- 
termediate time scales in a manner analogous to tempo- 
ral relaxations observed in glassy systems at low tem- 
peratures. Additionally, the time scale extracted from 
measurements of the stress relaxation curve C\ (t) is seen 
to diverge as l/vj. This time scale can be associated 
with the average lifetime of a highly-stressed/frequcntly- 



colliding chain in the system and thus the chains appear 
to become infinitely long-lived as v/ — » 0. 

It is worth noting that a similar search for a relevant 
length scale associated with these chains has been more 
difficult. Simple two-point spatial correlations of the col- 
lisional stress have not yielded any meaningful results and 
thus a higher-order correlation function such as the four- 
point correlations measured in supercooled liquids [T(| 
measured may be necessary. Initial calculations of spa- 
tial velocity correlations in our flowing system have indi- 
cated the existence of a length scale which is increasing 
as the flow velocity is decreased 19] . Preliminary mea- 
surements show that this length scale is of the order of 
the system size at the slowest flow velocities measured 
thus far. Therefore another important consideration is 
finite-size effects; while we have not explored this issue 
in this preliminary work, some finite-size scaling may be 
useful in accurately determining any measurable length 
scale pgj . 

As discussed above, the trap model picture may pro- 
vide a framework for describing the dynamics of our sim- 
ple granular flow system as well as other supercooled 
liquids. To explore this possibility the relevant quan- 
tity corresponding to the distribution of trapping times 
must be determined and measured, and the predictions 
of this type of coarse-grained dynamics quantitatively ex- 
amined. 
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